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ABSTRACT 

Cosmological gamma-ray bursts (GRBs) and their afterglows seem to result from 
dissipation of bulk energy in relativistic outflows, but their engine has not been 
unambiguously identified. The engine could be a young pulsar formed from accretion 
induced collapse with a dynamo amplified field. Elsewhere, we suggest that such a 
"Usov type" strong field pulsar may help explain the bimodal distribution in GRB 
durations. Here we discuss possible roles of a pulsar for the afterglow. We derive the 
expected bolometric luminosity decay. The extracted rotational energy could dissipate 
by shocks or by large amplitude electromagnetic waves (LAEMW). The simplest 
LAEMW approach predicts a slower decay in observed afterglow peak frequency and 
faster decay in flux than the simplest blast-wave model, though more complicated 
models of both can provide different dependences. LAEMW do not require the rapid 
magnetic field amplification demanded of the blast-wave approach because the emission 
originates from a nearly fixed radius. Different time dependent behavior of GRB and 
post-GRB emission is also predicted. Observational evidence for a pulsar in a GRB 
would make some GRB engine models, such as neutron star mergers and black holes 
unlikely. Therefore, the question of whether a pulsar is present is an important one 
even if it could drive a canonical fireball. 



Subject headings: gamma rays: bursts; pulsars: general; stars: magnetic fields 



1. Introduction 



Gamma-ray bursts (GRB) and their afterglows are likely the result of energy dissipation from 
a relativistically expanding outflow (e.g. Rees 1997a). Recent observations of GRB and their 
afterglows in optical and radio require a total GRB energy output ~ 10 52 A<1 ) erg where is the 
solid angle fraction for the opening of the outflow from which emission occurs (e.g. Waxman et al. 



- 2 - 



1997, Dar 1997, and references therein). The observation of an intervening absorber at redshift of 
0.835 (Metzger el. al 1997) in GRB 970508, suggests that this burst, if not all, are cosmological. 
Although cosmological fireball+blast-wave models provide somewhat successful fits to the GRB 
and afterglow emission (e.g. Wijers et al. 1997, Waxman 1997, Katz & Piran 1997), the engine 
driving the outflow is unknown. This gives GRB a unique place in modern astronomy as the only- 
known source whose engine does not yet point to a convincing paradigm. Even the energy budget 
is not fully constrained since the extent of beaming is unknown (e.g. Rhoads 1997). 

Magnetic fields are essential to GRB emission and afterglows in most viable models. Two 
paradigms for the GRB and/or the afterglow can be distinguished (e.g. Katz & Piran 1997). In 
the first, emission results from dissipation in an out-flowing fireball and magnetized blast-wave 
interactacting with the external medium, "external shock model." In the second, energy is 
transferred into particles well internal to the interface between the outflow and the ambient 
medium. Internal shock models are an example of the latter but the class is better described 
just by "internal" since shocks are not necessarily required. It appears that external models are 
not at work for the actual GRB, but may be fundamental for the afterglow. Whether internal 
processes may account for both is unknown. One example of a "non-shock" internal process 
is large amplitude electromagnetic waves (LAEMW), which can extract rotational energy (e.g. 
Michel 1984, Usov 1994) from strongly magnetized, millisecond (or sub-millisecond) pulsars (MSP) 
and can be applied to GRB. 

Several pulsar or compact object models (e.g. Usov 1992, Thompson 1994, Blackman et. al 
1996, Kluzniak & Ruderman 1997; Meszaros & Rees 1997, Yi h Blackman 1998) share the feature 
that power is first extracted in the form of Poynting flux and is later converted to gamma-rays and 
lower frequency photons. Since the Poynting flux energizes particles at large distances from the 
magnetized rotator, the required low baryon fraction (e.g. Fenimore et al. 1993) is maintained, 
and the gamma-ray transparency condition on the bulk Lorentz factor of the relativistic outflow 
(7 > 100) is ensured. The Poynting flux extraction of Meszaros & Rees (1997) comes from a 
rotating black hole+torus engine, so the extraction is transient, lasting until the torus falls into 
the hole (tens of seconds), not throughout the afterglow, like a pulsar. 

Pulsars formed from accretion induced collapse (AIC), with a dynamo amplified field, are 
appealing for GRB because two classes of such bursts can result, depending on whether or not the 
initial pulsar forms with spin above or below the gravitationally unstable limit (Usov 1992). We 
have suggested elsewhere (Blackman et al. 1996, Yi &: Blackman 1998, cf. Katz Sz Canel 1996) 
that these two classes might account for the bimodal duration distribution (e.g. Kouveliotou et 
al. 1993). The pulsar also provides an alternative to explain the GRB afterglows. Also, pulsar 
LAEMW emission mechanisms exist which do not require the amplification of magnetic field 
required in blast-wave shock models. 

We emphasize that even though the presence of a pulsar may be fundamental, the actual 
mechanism by which the pulsar ultimately transfers its energy into radiation is not well understood. 
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If a blast-wave is produced, the result may be very similar to transient rotator models (e.g. 
Meszaros & Rees 1997) formed from neutron star mergers (e.g. Narayan et al. 1992). But if 
observational evidence for a pulsar is found, GRB would probably not be the result of neutron 
star mergers since the compact object formed is likely too large to be a neutron star (Ruffert et 
al. 1997). The question of a pulsar presence is thus extremely important. 

In this paper we first summarize an AIC (cf. Yi & Blackman 1998; also see Dar et al. 1992, 
Thompson 1994) pulsar model. We then discuss how the electromagnetic energy might dissipate. 
We calculate features of an LAEMW afterglow and compare to the external shock fireball. We 
then conclude, suggesting features that could help indicate a pulsar presence. 



2. Pulsar Model and Bolometric Luminosity Evolution 

When a white dwarf reaches the critical Chandrasekhar mass ~ 1.4M Q through mass 
accretion, it collapses to a neutron star. The probable mass accretion rate leading to AIC is 
M^3x 10 18 <7/s (e.g. Livio & Truran 1992), and lasts for > 10 6 yr. (Such constraints disfavor 
binary systems where the secondary is also a white dwarf.) For a pre- AIC white dwarf with moment 
of inertia I w d = 10 51 gcm 2 and radius R w d = 10 9 cm, the magnetic field B W( [ and spin frequency 
£l w d are related to the post- AIC pulsar spin frequency and magnetic field by £l w d = £l*(I*/I w d) 
and B w d = B*(R*/R w d) 2 where i* = 10 45 gem 2 and i?* = 10 6 cm are the moment of inertia and the 
radius of the MSP. In order to create a MSP (fi* ~ lO 4 ,^ 1 ) with B* ~ 10 15 G by flux-freezing, the 
pre-collapse white dwarf must have $l w d ~ 10~ 2 s and B w d ~ 10 9 G. Such strong white dwarf fields 
have not been observed. Moreover, Yi & Blackman (1998) have shown that such AIC inducing 
magnetized accretion is not compatible with such a white dwarf. They find that the most likely 
AlC-produced MSP parameters satisfy ~ lO 4 ^ i^s -1 where B^n = B*/W n G. 

But field amplication is natural in the early phases of the neutron star, either linearly by 
differential rotation (Kluzniak & Ruderman 1997) or exponentially by dynamo amplification 
(Duncan Sz Thopmson 1992). The dynamo would undoubtedly be the more efficient mode if 
operating. In the Kluzniak & Ruderman (1997) approach (which could in principle be applied 
to any magnetized rotator) the growth of magnetic field from differential rotation leads to 
intermittent expulsion of field from the pulsar, possibly accounting for the spiky light curves of 
GRB. This requires linear magnetic field growth rather than (dynamo) exponential, otherwise the 
expulsions would be too rapid. 

Duncan & Thompson (1997) suggested that a young hot MSP is a favorable site for an auj 
dynamo due to vigorous convection driven by a large neutrino flux. The dynamo efficiency can 
be estimated by the Rossby number Nr = P*/t cotiv , where P* = 2^/0* is the MSP spin period 
and T conv is the convective overturn time scale at the base of the convection zone. If Nr ^ 1 
in a turbulent medium, an efficient dynamo can result. Duncan and Thompson (1992) give the 
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convective overturn time as T con v ~ 10 3 i ? 39 s where F39 is the convective neutrino heat flux in 
units of 10 39 erg/s/cm 2 . They show that a ~ 10 15 G large scale field can be produced. 

The total luminosity emitted from the young pulsar has two important terms. An 
electromagnetic (EM) term and a gravitational radiation (GR) term. When the initial pulsar 
spin exceeds a critical spin £l C rit, the GR luminosity term can dominate (e.g. Usov 1992). For a 
young and highly viscous pulsar, the GR drain does not have time to change the star's moment 
of inertia (Yi & Blackman 1998), so GR proceeds through spin-down, and the pulsar follows the 
track of an unstable Maclaurin spheroid (e.g. Chandrasekhar 1969). This means that the EM 
luminosity is depleted on the spin-down time scale, enabling the AlC-pulsar model to account for 
the bimodality of GRB durations (Blackman et al, 1996; Yi &: Blackman 1998, however, cf. Katz 
& Canel 1996). If the young NS initially followed the Jacobi rather than Maclaurin track, it would 
instead spin up, and would not naturally lead to a bimodal distribution. However, the GR from 
these two paths may be measurably distinguishable (Lai & Shapiro 1995). 

Assuming £l C rit > ^dynamo, where Qdynamo is the spin required for dynamo action, there 
are two types of pulsars which can form GRB: [1] Supercritical strong field rotator (SPS) with 
SI* > Qcrit > ^dynamo, and [2] Subcritical strong field rotators (SBS) Q cr it > ^* > ^dynamo- Since 
these classes have similar rotational frequencies, they likely originate from similar pre-collapse 
conditions so the numbers of the two classes should be similar. The EM luminosity available for 
particles (quadrupole may also be important) is 

L obs = fL ern ~ 10 51 (//0.1) (A^/O.Ol)- 1 (O./10 4 s' 1 ) 4 {R* /10 6 cm) 6 (£„ /10 15 G) 2 erg/s, (2-1) 

where / is the efficiency. The luminosity from GR is given by 

L gw = 32Ge 2 / 2 f} 6 /5c 5 = 7 x lO 55 (M*/M )(i?*/lO 6 cm) 4 (O*/lO 4 s~ l fe 2 erg/s, (2-2) 

where e is the eccentricity. Both SPS and SBS emit EM radiation. However for SPS, the 
spin-down, and hence the luminosity e-folding decrease, occurs on a time scale r gw ^ Is whereas 
in SBS, the spin-down and luminosity decrease occur on a time scale r em ~ 10 2 s. SPS may be 
naturally related to the short GRB bursts and SBS to long bursts of the observed (Kouveliotou et 
al. 1993) bimodal duration distribution. 

Using (2-1) and (2-2), we can solve for the time evolution of the PSR spin for L gw and L em 
dominated regimes respectively. Setting the kinetic energy loss L kin = (l/2)M*Rln*(dn*/dt) 
equal to L gw , we use (2-2) to solve for fi*. Separating variables and integrating gives 

fi, = 10 4 (fVl0 4 s _1 )[l + 2.8 x 10 3 (^ /10 4 s- 1 ) 4 (M*/M )( J R/lO 6 cm) 2 e 2 ((t - t )/ls)}- 1/4 s~\ 

(2-3) 

where Qo = f2*(t = to) is the spin rate of the initially formed pulsar at time t = to. At t = t eq , 
the spin evolution for SPS bursts changes from being GR dominated to EM radiation dominated. 
This time is likely of order r gw since e decreases rapidly after that, and L gw becomes unimportant. 
For the SBS, the spin down is always dominated by Lem- 
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For the SBS bursts and for SPS bursts after t eq , we set (2-1) equal to L^in and integrate, so 

ft, = 10 4 (1V10 4 s- 1 )[l + 2x 10- 2 (O /10 4 s- 1 ) 2 {R/lQ) % cm)' i {B^/l^G) 2 {{t-t ()>em )/ls)}- 1 ' 2 a" 1 . 

(2-4) 

where to,em equals to for the SBS or t eq for the SPS. 

By plugging (2-3) and (2-4) back into (2-1) we can find the observed bolometric luminosity 
time dependence for GRB and afterglow. For the GR dominated phase 

L obs ~ 10 51 (//0.1)(A*/0.01)- 1 (12*/10 6 cm) 6 (B,/10 15 G) 2 (f2o/10 4 s" 1 ) 4 (2-5) 

x [1 + 2.8 x 10 3 (O /10 4 s' 1 ) 4 (M„/M Q )(R/W 6 cm) 2 e 2 {{t - ^/ls)]- 1 erg/s, 
while for the EM dominated phase 

Lobs ~ 10 51 (//0.1)(A*/0.01)- 1 (i2./10 6 cm) 6 (B,/10 15 G) 2 (f2 0l em/10 4 s' 1 ) 4 (2-6) 

x[l + 2 x 10- 2 (O , em /10 4 s- 1 ) 2 (/2/10 6 cm) 4 (B»/10 15 G) 2 ((t-*o 1 em)/la)]" 2 erg/s, 
where fio,em is the spin at the beginning of EM domination at t = to >em . 

Note that in the very early phase of the SBS bursts, the spin-down driven emission remains 
approximately constant and does not begin to decay as a power law until the critical time t ~ 50s 
from (2-6). Whereas for the SPS bursts, from (2-5), a t~ l decay occurs immediately, followed by a 
t~ 2 decay from the L ern dominated phase of (2-6). Thus SBS are in qualitative accordance with 
the data analysis of Fenimore (1997) in the sense that during the GRB the burst decays more 
slowly than in the subsequent afterglow. This is just an example of how the same internal process 
can account for different time evolutions of the burst and afterglow. 



3. Peak Frequency and Associated Luminosity Evolution from LAEMW 

The magnetic field in the pulsar engine acts as a drive belt, extracting energy as Poynting flux 
to a large distance where it is then converted into particles. (The magnetic field plays a similar 
role in the Meszaros & Rees (1997) model.) The extracted energy could subsequently dissipate in 
large amplitude electromagnetic waves or within a fireball+shock. How this energy is dissipated 
or converted into particles is not clear and may depend on the proton fraction of the outflow. If 
there is a negligible proton fraction, then LAEMW can propagate and accelerate pair plasma blobs 
(Asseo et al. 1978, Blackman et al. 1996). We discuss this and the associated afterglow below. 

Avoiding the runaway pair production that would make a blob optically thick to gamma rays, 
requires 7 > 100 (e.g. Krolik & Pier 1991). Such pair plasma blobs moving along the magnetic 
dipole axis with 7 > 10 4 might be produced in PSR magnetospheres (Usov 1994) by the LAEMW. 
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The LAEMW propagate outside a radius, r//, where the density drops below that which can 
sustain a Goldreich-Julian (Goldreich & Julian 1969) charge density, and thus where flux-freezing 
and force-free conditions are broken. This gives (Usov 1994) 

r ff ~ 7.5 x 10 13 (//0.1) 1 / 4 (A*/0.01)- 1 / 4 (B,/10 15 G) 1 / 2 (n # /10 4 s- 1 ) 1 / 2 cm, (3-1) 

where we have included the beaming. The associated pair plasma number density is 

n ff ~ 1.5 x 10 8 (//0.1)- 3 / 4 (A$/0.01) 3 / 4 ( J R*/10 6 cm)(^/10 15 G) 1 / 2 (^/10 4 s- 1 ) 5 / 2 cm~ 3 . (3-2) 

Electron acceleration at rff is characterized by aff, defined by (Usov 1994; Michel 1984) 

Off = L dip /{mc 2 Nf f ) ~ 5.3xl0 7 (//0.1)- 3/4 (A$/0.01) 3/4 ( J R,/10 6 cm) 5 ( J B,/10 15 G) 1/2 (O*/10 4 s- 1 ) 1/2 , 

(3-3) 

where Nff = ^irr 2 ^criff is the electron flux and m e is the electron mass. By solving the equations 
of motion for a particle in a pulsar wind zone subject to electromagnetic forces, it has been shown 
that relativistic electromagnetic waves of frequency 0* can accelerate pair plasma to (Michel 1984, 
Asseo et al. 1978) 7 ~ ajf ~ 10 6 (i?*/10 6 cm) 10 / 3 (S*/10 15 G) 1 / 3 (^*/10 4 s" 1 ) 1 / 3 , and the resulting 
emission is beamed within solid angle ~ 7~ 2 from the direction of wave propagation (Asseo et al. 
1978). The characteristic emitted frequency of the synchro-Compton radiation is (Michel 1984) 

v c ~ ln{rf f /R,)n,a 2 ff ~ 5xl0 20 (//0.1)- 3 / 2 (A$/0.01) 3 / 2 (^,/10 4 s" 1 ) 2 ( J R*/10 6 cm) 10 ( J B,/10 15 G) Hz. 

(3-4) 

Using (3-4), (2-3) and (2-4), and assuming only fi* changes with time, we obtain different 
dependences of v c {t) for the GR and electromagnetic dominated phases. For the GR phase 

u c = zy cO [l + 2.8xlO 3 (^ cO /5xlO 2O ^)(i? ; ,/lO 6 cm)- 18 (^/lO 15 ) 2 (M/M )e 2 ((t-to)/ls)]" 1/2 , (3-5) 

where v c o = v c {Q* = &o), while for the electromagnetic dominated regime 

vc = v c0 ,e m [l + 2 x 10- 2 (^ c0 , em /5 x W 20 Hz)(R*/10 6 cm)- g (B*/10 15 )({t - t^ em )/ls)\-\ (3-6) 

where is c o,em = ^ c (^* = ^o,em)- Note that when (3-5) applies, the total emission is actually 
dominated by GR which peaks at 10-100 Hz (Lai & Shapiro 1995). 

We can also compute the time decay of the luminosity (or flux) at the peak electromagnetic 
emission frequency v c . Assuming the form L u = L Uc (v /v c )~ a , with a > 1, we have L b s ~ L Uc v c . 
Using L Vc ~ u~ 1 L f )S , and for the GR dominated phase using (3-4) and (3-5) in (2-1) we get 

L Vc ~ 10 51 (//0.1)(A$/0.01)- 1 (^/10 6 cm) 6 ( J B,/10 15 G) 2 (O /10 4 s" 1 ) 4 (3-7) 

x [1 + 2.8 x 10 3 (S!o/10 4 s- 1 ) 4 (M*/M Q )( J R/10 6 cm) 2 e 2 ((t-t )/ls)]~ 1/2 erg/s. 
For the electromagnetic dominated phase using (3-4) and (3-6) in (2-1) we get 

L Vc ~ 10 51 (//0.1)(Ad»/0.01)- 1 (^/10 6 cm) 6 (^/10 15 G) 2 (^o, em /10 4 s" 1 ) 4 (3-8) 
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x[l + 2 x 10- 2 (O ,em/10 4 s- 1 ) 2 (R/ 10 6 cm) 4 (B*/W 15 G) 2 ((t-t , em )/ Is)}' 1 erg/s. 

The evolution of the critical frequency and associated luminosity is less certain than the evolution 
of the bolometric luminosity when the latter is directly related to pulsar spin-down. Nevertheless 
we can still speculate that the proton fraction in the outflow may be essential in determining 
what kind of dissipation occurs and the LAEMW case discussed above is most relevant when the 
outflow contains very few. Simulations of Hoshino et al. (1992) show that collisionless relativistic 
shocks produce a significant non-thermal lepton tail only when the protons comprise > 5% of the 
material. In the proton enriched case, shock acceleration of electrons might be the most likely 
GRB producing mechanism. The simplest fireball evolution (e.g. Waxman 1997) leads to a t~ 2 
bolometric luminosity dependence, just like the LAEMW, but predicts v c oc i -3//2 and F Vc oc i -1 / 2 . 
This does seem to be consistent with several afterglows (Wijers et al 1997; Waxman 1997). 

The characteristic "smoothness" of the outflow Lorentz factor V may also be important in 
distinguishing what kind of dissipation occurs. For ejecta emitted at two different times with two 
different V satisfying Y2 > T\ » 1, internal GRB producing shocks could form at an observer 
frame distance of c|*2 — iilTi^ (Rees & Meszaros 1994) and may dissipate the energy before other 
internal mechanisms like LAEMW acceleration occur. 



4. Discussion 

A generic afterglow may not clearly reveal the GRB engine, much like supernova remnant 
physics is disconnected from the physics of stellar collapse. A pulsar GRB engine could lead to a 
fireball or alternative acceleration mechanisms like LAEMW that dissipate a large fraction of EM 
energy before shocks form at a nearly fixed radius. One advantage of this LAEMW approach is 
that there is no need for the extremely rapid magnetic field amplification in the expanding outflow 
required of the fireball for both external and the internal shock models (e.g. Rees & Meszaros 
1994): For LAEMW, the radius rjj and density at which the particle acceleration and radiation 
can occur remains nearly constant while the characteristic emission frequency decreases from its 
sensitivity to pulsar spin down. Alternatively, the pulsar could also provide some of the required 
magnetic field, if its role were instead to generate a blast-wave. 

The simplest LAEMW picture predicts slower afterglow critical frequency decays than a 
blast-wave. However, variations in the blast-wave models might lead to slower decays as it is not 
necessarily the case that most of the outflow energy is carried by the flow moving with the largest 
Lorentz factor (Rees & Mezaros 1997). To complicate things further, both blast-waves and LAEM 
could result from pulsars. The proton fraction of the outflow may be important in determining 
the primary mode of particle acceleration and emission. In short, variations of energy dissipation 
mechanisms in strongly magnetized pulsar outflows could lead to a variety of critical frequency 
and peak flux evolutions. This makes it hard to distinguish between a pulsar engine vs. that of a 
transient Poynting flux extraction model. 
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However, whether a pulsar lies in GRB engines is important because the presence would 
be incompatible with the resulting compact object mass from neutron star merger models and 
neutron star equations of state (Ruffert et al. 1997). Evidence for pulsars in GRBs would 
threaten the merger paradigm. Several possibilities that could hint at a pulsar are: [1] The 
pulsar rotation itself may be a source of ~ 10 _3 s. variability during the burst. For SBS (the 
longer) bursts, the characteristic time scale of variability after 50s or so from (2-6) goes as 
t v oc fi^r 1 oc i 1 / 2 . This trend may be measurable in X-rays. For SPS (the shorter) bursts, 
the trend is t v oc oc i 1 / 4 , immediately, and then t v oc Q^ 1 oc t 1 / 2 after r gw . [2] Future 
evidence in the radio band (Rees 1997b) like violent fluctuations or coherent emission (e.g. 
Melrose 1993). For GRB970508, the predicted bolometric luminosity from the t~ 2 decay after 
t ~ 50days is about equal to that required of the radio afterglow luminosity for GRB970508 which 
is L obs ~ W il {v/5GHz)(F u /W- 26 erg/s/cm 2 )(d/5Gpc) 2 erg/s (Waxman et al. 1997): Using (2-6) 
for SBS, we have L obs ~ 10 41 erg/s at t = 50 days. However, the decreasing scintillation amplitude 
suggests radio source expansion (Goodman 1997, Waxman et al. 1997), so here any pulsar's role 
could have been mainly to drive a blast-wave, rather than LAEMW. [3] The SPS class of GRBs 
would be accompanied by GR absent from the longer SBS class. The detection of a significant GR 
excess from short GRBs would give support to an AIC pulsar engine model. The GR signature for 
various evolutionary tracks of secularly unstable AIC pulsars are unique and possibly measurable 
(Lai & Shapiro 1995). [4] The SPS bursts should show a t~ l decay in bolometric luminosity 
followed by a t~ 2 decay after the gamma-ray phase. The SBS bursts should show almost no decay 
during the bursts and then a t~ 2 decay after the GRB phase ends. This follows from (3-7) and 
(3-8). [5] Observations might reveal a faint residual pulsar or X-ray accreting source well after the 
optical afterglow fades. An AXAF view of an afterglow could be revealing. 

Finally, the AIC rate for pulsar formation is unknown, but could be comparable to the local 
supernova rate, i.e. 10 _2 yr _1 per galaxy, the observed GRB rate of ~ 10 -6 — 10~ 5 yr~ 1 per 
galaxy can be amply explained by the AIC model with modest beaming. These numbers are not 
inconsistent with the possibility that GRB signature all AIC events. 
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